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We have calculated analytically the mean value and the variance of the number of bonds on the 
lattices of dimension d for the given occupation of sites. We consider both kinds of site occupation: 
with the fixed concentration n s of occupied sites and with the probability p for a site to be occupied. 
We found that the variance are smaller in the former case and not depends from the dimensionality 
of the hypercubic lattice. Whereas in the last case it grows with the space dimensionality. The ratio 
of variances are quite different in the limit of p — * 1. Finally, we demonstrate the relevance of the 
level of fluctuations on the precision of energy calculations for the Ising model in the Monte Carlo 
simulations. 



I. INTRODUCTION 

There are two ways to generate on the lattice a sample with the randomly occupied sites H . The first one, with 
the fixed number s of occupied sites. This is the sampling from a 'canonical ensemble' for which in this paper we will 
use the name s-sampling. The second way, with a fixed probability p for a given site to be occupied, i.e. sampling 
from a 'grand canonical ensemble', or p-sampling. It is known that both methods lead in the 'thermodynamic' limit 
of infinite lattice size to the same results for the free energy per Ising spin being placed onto occupied sites of the 
lattice J|. 

If one choose the probability p for the site to be occupied the concentration of bonds between the occupied sites 
on the infinite lattice are nss — p 2 ■ We named this bonds as S'S'-bonds. Similarly, the concentration of bonds 
between non-occupied sites (iViV-bonds) is rijv./v = (1 — p) 2 and the concentration of bonds between occupied and 
non-occupied sites (SiV-bonds) are usn = 2p(l — p). Clearly, the sum of the three concentrations nss + n<NN + nsN 
is the probability that each bond belongs to one of the three types and equal to unity. 

In this paper we asking the following questions: i) how fluctuates the concentration of bonds on the hypercubic 
lattices depending on the concentration of occupied sites? ii) how this fluctuations depend on the way of sampling? 
iii) how this fluctuations depend on the lattice size? 

We found analytically the answers on the all of the stated questions. 

We start our analyses with the simple case of one-dimensional periodic lattice in Section ||. In Section ||| we 
calculated the mean value and the variance of the number of bonds for the square lattice using the results of Section [n|. 



In the following Section IV we genelarize our results for the d-dimensional hypercubic lattice. We conclude our paper 
with the discussion of results in Section 

II. ONE-DIMENSIONAL LATTICE 
A. s-sampling 

Let us consider an one-dimensional lattice with L sites and periodic boundary (i.e. ring). 

We are interested in the probability P(b; s, L) to have b bonds between s occupied sites. The number of occupied 
sites s, the number of bonds between occupied sites b and the number of clusters of occupied sites k are connected 
by relation 



(1) 



Then, the probability could be easily written as 



LC k s zlC 



P(b;s,L)= fc C r , ( 2 ) 
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where combinatorial coefficients C^Zi is the number of dispositions of s sites among the k clusters, C£Zs_i - the same 
for L — s empty sites, C S L - number of configurations with the s occupied sites from L possibilities, and, finally, L/k - 
the number of positions of k clusters on L sites. The probability are normalized 

s— 1 s— 1 j /~<s— 6— Isis— 6— 1 

Ep(m,d=E (/-6)^- l =L ( 3 ) 

6=0 6=0 V ' L 

Now, the mean value of bonds (b) are 

(b)=Y / bP(b; S ,L)='^-^- (4) 

6=0 

and the mean density of bonds are given by L » 1 

n b =(b)/L^n 2 s -n s (l-n s )/L, (5) 

where n s = s/L is the density of occupied sites. 

The variation of the number of bonds Ab — (6 2 ) — (b) could be calculated in the same manner 

and the variance of the density of bonds Ani, are equal to 

An b k n^(l - « s ) 2 - - n s )(4yi2 - 4n s + 1)/L (7) 

for L >> 1. 

B. p-sampling 

The probability distributions of the some quantity A in the case of s-sampling P(A; s) and those of the case of 
p-sampling P(A;p) are connected by 

L 

P(A;p)^P(A; S )C£p s (l-p) L - s , (8) 

s=0 

therefore the probability to obtain exactly b bonds between two sites occupied with the (independent) probability p 
on the lattice of size L are 

P(b;p, L) = J2 P) L ^ ( 9 ) 

s=0 \ L J 

with fc given by eq. (Q) and the mean value of the bonds (&)' in the case of p-sampling will be given by 

s-l L 

(b)' = J2bP(b;p,L)^J2( c lP s ( 1 -P^ s ( b ))= L P 2 > ( 10 ) 

6=0 s=0 

(where prime sign ' denotes the p-sampling averages) and the corresponding variation Ab' are given by 

L 

Ab' = J2(CIP S (1-P) L - S (A6(6) 2 )) -(b)' 2 = V(l-p)(l + 3p) (11) 

s=0 

for L > 4. 

We have to stress here, that the mean density of the number of bonds n p — p 2 and their variance An p = p 2 (l — 
p)(l + 3p) do not contain any finite size dependences in contrast with the case of s-sampling behaviour in (||) and |?|). 
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III. SQUARE LATTICE 



Let us calculate the mean value of the number of bonds and it's variation for the square lattice with periodic 
boundary conditions. 



A. p-sampling 

On the square lattice with periodic boundaries, there are L 2 vertical and L 2 horizontal bonds. The number of 
bonds in each of L columns and each of L rows are the random variables. Their mean values (b)^ are given by ( |l(i| ) 
and their variances by (pd|). Thus, the mean value of the number of bonds on 2d lattice are 

2L 

<^ 2 = E^ = 2 ^ 2 - ( 12 ) 

i=l 

and their variance are given by |3) 

2L 2L 

Ab' d=2 = J2 Ab ' + 2 H Cov (< 6 >'< . ( 13 ) 



1=1 l,j 

where i < j. 

It is clear that the column-column (row-row) covariances are equal to zero, and the column-row covariances are the 
same for all columns and rows, and 

Ab' d=2 = 2LV(1 - p)(l + 3p) + 2L 2 Cov {{b)' column , (&>' ro J . (14) 
We have to compute now the column- row covariance Cov ((b)' column , {b)' row ) 

Cov {(Kolumn > ( b )'ro W ) = Wcolumn ( b )'row) ~ ( b )' 2 ■ ( 15 ) 

The first average are splitted into two parts. The first one, corresponds to all cases of distribution of occupied sites 
on columns and rows with the constraint that the intersection of row and column a;* are occupied. The second one, 
corresponds to the states with the unoccupied intersection site. 

(( b )'coiumn ( 6 )r<™) = P ( b \ x * = occupied) 2 + (1 - p) (b\x* = unoccupied) 2 (16) 

We omit here the algebra leading to the final expression for the variation of the number of sites on the two- 
dimensional lattice with the given probability p for the site to be occupied 

Ab' d=2 = 2L 2 p 2 (l-p)(l + 7p). (17) 

Details will be published elsewhere B. 



B. s-sampling 

We calculate also the mean value of the number of bonds on the square lattice with the linear size L and with the 
fixed number s of occupied sites 



<^ 2 = ^^ 



and the variance are 



s(s - l)(L 2 - s)(L 2 - s - 1)(2L 2 - 10) 



Abd = 2 ~ (L 2 -l) 2 (L 2 -2)(L 2 -3) • (19) 
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IV. HYPERCUBIC LATTICES 



The mean value of the number of bonds on the d-dimensional hypercubic lattice with the probability p for the site 
to be occupied are 

(b)' d = dL d p 2 . (20) 

The varian ce could be calculated in the same manner as in the two dimensional d = 2 case described above in the 
Section IlII A 



Ab' d = dL d p 2 (l - p){\ + {Ad - l)p) . 



(21) 



The explicit mean value and the variance of the number of bonds for the case of the fixed number s of occupied 
sites could be calculated using ( p0| ) and ( |2l] ) with the proposition that the mean value are quadratic in s. The result 
are 



ds(s - 1) 
L d - 1 



(22) 



and 



A6 d = 



ds(s - l){L d - s){L d -s - \){L d -2d -I) 
(L d -l) 2 (L d -2){L d -3) 



(23) 



for L > 2. 

The expressions (pfl-refanydbs) we checked by direct enumeration in d — 2 for L = 4, 5 and in d = 3 for L = 3 



V. DISCUSSION 



We do not find any evidence for the absense of self-averaging [gj^] in the number of bonds. In the case of p-sampling 
the normalized variance 

fit , = A| = (l-»d (24) 

d (by* p 2 L d 

inverse proportional to the lattice volume. Therefore the number of bonds are the strong averaging quantity 
The normalized variance in the case of s-sampling 



Ab d (l~n s ) 2 
(bf d ^ dL- d n s (n s - l/L d ) 



are nonmonotonic in L and the strong self-averaging take place for L > (l/n s ) 1 ^ d . 
The relative variance for two cases 

Ab[i = P 2 (l + (4d-l)p)(l-p) 
A6d (n s (l-n s )f 

shows that in the limit of p — * 1 and n s -> 1 the ratio diverges. 

It is possible to give the physical interpretation to the last expression. Let us place Ising spin at each occupied site. 
Then, the energy Eq at zero temperature will be equal to the number of bonds betwen the sites occupied by spins. 
The magnetization M are equal to the number of occupied sites. Then, the ratio ( |26| ) are equal to the ratio N p /N s 
of the number of samples (N p for the p-sampling and N s for the s-smapling) over which we have to average energy 
Eq in order to get the same dispersion uEq = -J Ab/ (N x — 1) of the energy Eq, where N x are either N p or N s . 

Finally, we could propose the third way to generate the samples, sb-sampling: get random samples with the fixed 
number of sites s and choose among them only those which contain the "right" number of bonds b = ds 2 . The 
results of Monte Carlo simulations for site diluted Ising model on square lattice had show Q that the dispersions 
of termodynamic quantaties like energy, specific heat and magnetic susceptibility are quite smaller for the proposed 
sb-sampling in comparison with both, s-sampling and p-sampling at all temperatures. 



4 



VI. ACKNOWLEDGMENTS 



Authors are thankfull to K. Binder for the kind discussion and suggestion to write this manuscript. The work is 
partially supported by grants from RFBR (99-02-18412), NWO and INTAS. OA.V. is grateful to Landau stipendium 
committee (Forshungzcntrum/KFA, Julich) for support. 



[1] L.N. Shchur, Incipient Spanning Clusters in Square and Cubic Percolation, in Springer Proceedings in Physics, "Computer 
Simulation Studies in Condensed Matter Physics XII", Eds. D.P. Landau, S.P. Lewis, and H.B. Schh'uttler, (Springer- 
Verlag, Berlin, 2000) 

[2] R.B. Griffiths and J.L. Lebowitz, J. Math. Phys. 9 (1968) 1284 

[3] W. Feller, An introduction to probability theory and its applications (John Wiley & Sons, New York, 1970) 
[4] O.A. Vasilyev, to be published. 

[5] B. Derrida and H. Hilhorst, J. Phys. C: Solid State Phys. 14 (1981) L539 

[6] K. Binder, D.W. Heermann, Monte Carlo Simulations in Statistical Physics (Springer-Verlag, Berlin, 1997) 
[7] S. Wiseman and E. Domany, Phys. Rev. E 52 (1995) 3469; Phys. Rev. Lett. 81 (1998) 22; Phys. Rev. E 58 (1998) 2938 
[8] A. Aharony and A. B. Harris, Phys. Rev. Lett. 77 (1996) 3700; A. Aharony, A. B. Harris and S. Wiseman, Phys. Rev. Lett. 
81 (1998) 252 



5 



